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ABSTRACT 

The effects of anisoplanatism on the adaptive optics point spread function are in- 
vestigated. A model is derived that combines observations of the guide star with an 
analytic formulation of anisoplanatism to generate predictions for the adaptive optics 
point spread function at arbitrary locations within the field of view. The analytic for- 
mulation captures the dependencies of anisoplanatism on aperture diameter, observing 
wavelength, angular ofi'set, zenith angle and turbulence profile. The predictions of this 
model are compared to narrowband 2.12 /im and 1.65 //m images of a 21 arcsec binary 
(m^=7.3, 7.6) acquired with the Palomar Adaptive Optics System on the Hale 5 me- 
ter telescope. Contemporaneous measurements of the turbulence profile made with a 
DIMM/MASS unit are used together with images of the primary to predict the point 
spread function of the binary companion. Predicted companion Strehl ratios are shown 
to match measurements to within a few percent, whereas predictions based on the iso- 
planatic angle approximation are highly discrepant. The predicted companion point 
spread functions are shown to agree with observations to 10%. These predictions are 
used to measure the diff'erential photometry between binary members to an accuracy 
of 1 part in 10^, and the differential astrometry to an accuracy of 1 mas. Errors in 
the differential astrometry are shown to be dominated by differential atmospheric tilt 
jitter. These results are compared to other techniques that have been employed for 
photometry, astrometry, and high contrast imaging. 

Subject headings: instrumentation:adaptive optics, techniques: high angular resolution, 
techniques: photometric 



1. Introduction 

An adaptive optics system senses phase aberrations arising from atmospheric turbulence using 
observations of a guide star, and compensates these aberrations by applying a correction to an 
adaptive mirror. This compensation is valid in the direction of the guide star, but degrades with 
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angular offset from the guide star due to anisoplanatism. This effect arises due to the shearing 
between the columns of turbulent atmosphere traversed by light from the guide star and light from 
a target at a finite angular offset. 

The phase aberrations that arise from anisoplanatism depend on a number of parameters. The 
errors grows with angular offset from the reference source, so that image quality and Strehl ratio 
degrade with increasing angular offset. The vertical distribution of turbulence has a strong effect 
on the degree of error, with higher altitude turbulence generating larger errors due to the larger 
geometrical shear. This dependence varies in time as the distribution of atmospheric turbulence 
over the telescope evolves. The anisoplanatic error grows with zenith angle, since one sees more 
turbulence along line of sight to the object. The error is also a strong function of aperture diameter 
and observing wavelength. 

This large number of dependencies yields a rich phenomenology for the adaptive optics point 
spread function (PSF). The parameter space is so large that it is not practical to integrate long 
enough to attain the stochastically averaged PSF. The variability in the adaptive optics PSF is a 
serious impediment to the quantitative interpretation of observations, and can limit the precision 
in astronomical applications involving photometry, astrometry, crowded field imaging, and high 
dynamic range imaging of extended objects. Integral field unit spectroscopy of extended objects is 
another application in which the dynamic range of observations may be seriously compromised by 
PSF variability. 

Multiconjugate (Beckers 1988) and multiobjcct (Hammer et al. 2004; Ellerbroek et al. 2005) 
adaptive optics architectures have been proposed that aim to directly overcome the effects of aniso- 
planatism. These architectures use multiple guide stars distributed over a finite field of view and 
employ tomographic algorithms to estimate the three dimensional volume of atmospheric turbu- 
lence. The algorithms rely on a knowledge of the angular offsets among the guide stars, and use 
the statistical correlations induced by these angular offsets in order to effect the turbulence esti- 
mation. In this sense, these algorithms employ anisoplanatism in order to overcome the effects of 
anisoplanatism. Anisoplanatism in a single conjugate adaptive optics system forms the limiting 
case of this more complex problem. On-sky tests with existing adaptive optics systems can serve 
an important role in validating our understanding of anisoplanatism and increasing our confidence 
in the success of tomography. 

Significant efforts have been directed towards overcoming the effects of anisoplanatism in ex- 
isting adaptive optics systems. One line of investigation has aimed to extract estimates of the PSF 
from observed data. For a small target field, the PSF may be assumed to be field independent. 
Observations of such fields that contain multiple point sources may be deconvolved using a refer- 
ence PSF selected from the data or by solving for the optimal PSF during deconvolution. These 
approaches have been applied to crowded stellar fields (Diolaiti et al. 2000; Christou et al. 2004) 
and planetary objects (de Pater et al. 2004). In an effort to calibrate wider fields, Steinbring et 
al. (2002) measured the effects of anisoplanatism on the PSF using observations of crowded fields. 
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and then applied these measurements to other fields of interest. This type of technique is unable 
to capture any temporal evolution of the turbulence profile that occurs between the observations. 
Another approach has employed a parameterized model of the effects of anisoplanatism on the 
adaptive optics PSF, extracting these parameters during the deconvolution procedure (Flicker k. 
Rigaut 2005). 

The above techniques did not attempt to use any independent information about the tur- 
bulence profile in estimating the adaptive optics PSF. Other researchers have incorporated such 
measurements into the PSF estimates. Voitsekhovich et al. (1998) computed the structure function 
due to residual phase aberrations from anisoplanatism, and used this to evaluate the Strehl ratio 
as a function of angular offset from the guide star for different turbulence profiles and aperture 
diameters. Fusco et al. (2000) derived an expression for the adaptive optics optical transfer func- 
tion that captured the depenedencies of anisoplanatism on the above parameters. These authors 
carried out 850 nm observations of two binaries with the ONERA adaptive optics system and used 
a turbulence profile measured from a balloon flight together with the analytical expression for the 
OTF to predict the binary companion PSF. Weifi et al. (2002) performed K band adaptive optics 
observations of a binary using the ALFA adaptive optics system while at the same time measur- 
ing the turbulence profile using a scidar instrument. These authors compared the Strehl ratio 
degradation due to anisoplanatism measured from the binary image data with that expected from 
the measured turbulence profiles. These experiments have shown a promising level of agreement 
between predictions and measurements. 

Recent work in the area of turbulence monitoring has yielded an automated set of equipment 
capable of delivering real time estimates of the turbulence profile on minute timescales (Tokovinin 
et al. 2005). This equipment is based on the combination of a Multi- Aperture Scintillation Sensor 

(MASS) (Kornilov et al. 2003) and a Differential Image Motion Monitor (DIMM) (Vcrnin & Munoz- 
Tunon 1995). This turbulence monitoring equipment is now being used at a number of different 
sites, and has been employed by the Thirty Meter Telescope project for its site testing program 
(Skidmore et al. 2004). As part of this program, a set of this equipment has been installed at 
Palomar Observatory. 

This paper describes an experiment in which measurements of the turbulence profile from 
the DIMM/MASS equipment and short exposure images of a 21 arcsec binary from the Palomar 
Adaptive Optics System on the Hale 5 meter telescope were acquired contemporaneously over the 

course of several hours. Section 2 presents the analysis of anisoplanatism used for this experiment. 
Section 3 describes the observations, while Section 4 and 5 describe the analysis and results of the 
experiment. 
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2. A Model for the Field Dependent AO PSF 

2.1. Model Formulation 

Consider an adaptive optics system that measures the wavefront phase aberrations 4>a{r) in 
the direction of a guide star and compensates these aberrations using an adaptive mirror. Here f is 
a vector in the pupil plane of the telescope. An adaptive optics system is not capable of effecting a 
perfect correction, and residual phase aberrations (pa (^) will be present in the guide star wavefront 
after adaptive compensation. The wavefront phase aberrations 06 (r) in another direction on the sky 
differ from those of the guide star due to anisoplanatism. In this direction, the residual aberrations 
that remain after compensation by the adaptive optics system are 

Acf>{r} = </,,(f)-0„(r) + 4(r) (1) 

This simple model does not account for time delay in the application of the adaptive correction. 
During this delay the wind carries turbulence past the telescope aperture, causing a servo error in 
the guide star wavefront. This turbulence evolution also induces correlations between the adaptive 
correction and the wavefront aberrations in other directions on the sky (Tyler 1983). Under cer- 
tain observing conditions, this effect can lead to a situation where the Strehl ratio can be higher 
downwind of the guide star (Lena & Lai 1999). In this analysis, this wind induced anisoplanatic 
effect is neglected. 

The structure function for a random process ip{r) is 

D^{ri,r2) = ({V(ri)-V'(r2)f) (2) 
Using Equation 1, one can write the structure function for the residual phase A(p (r) as 

Da<i> (n , f2 ) = Dapi (r 1 , f2 ) + D^^ (ri , r2 ) + (3) 

2 <^|0apl(rL)0a(ri) + ^Z!>apl(r2)0a(^2) - 0apl (n ) 0a (^^2 ) - 0apl(^2)0 a iri)}) 

where 0api(^) = 0fc (^) — 4>a (^ is the component of the residual phase arising from anisoplanatism. 
The four cross terms in Equation 3 represent statistical correlations between the residual phase 
errors from the adaptive optics correction and those from anisoplanatism. A number of error terms 
that arise in adaptive optics systems have no correlation with anisoplanatic errors. Examples 
include measurement errors, aberrations in the optical system, and non-common path errors. In 
contrast, residual fitting error and servo error are weakly correlated with anisoplanatic errors. While 
the relative importance of these cross terms will depend on the error budget of the AO system, 
in many circumstances these terms will be small relative to the structure functions in Equation 3. 
Here, these cross terms are neglected. 

With this approximation, the optical transfer function may be written as (Goodman 1985) 

OTF{r) = J dsexp!^~[D^pi{s,r + s) + D^^{s,f+s)\^W{s)W{r + ^ (4) 
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where W{r) is the pupil function. An issue of particular importance in evaluating the OTF is 

whether the structure functions of the two random processes are stationary over the pupil plane. 
For such a process, D^{fi,f2) = D,p{r2 — n)- As will be shown in the next section, -Dapi(ri,r2) is 
in fact stationary. This permits the OTF to be written as 

OTF{r) = exp|-iz)api(rl|y"dsexp|-^Z^^Js,r + sl|w(slW^(r + s-) (5) 

Note that the OTF has factored into a term that describes the anisoplanatic errors and a term that 
describes the residual errors in the direction of the guide star. 

A relationship similar to Equation 5 was derived by Fusco et al. (2000) , and represents a very 
important result. One may formulate the OTF in any direction on the sky from the product of the 
guide star OTF and an anisoplanatic transfer function formed from -Dapi(^)- This factorization is 
of considerable use in the formulation of a practical scheme for evaluating the OTF, as discussed 
below. 



2.2. The Anisoplanatic Structure Function 



Prediction of the OTF in Equation 5 requires an evaluation of the anisoplanatic structure 
function Dg^pi{r). This function may be computed using a semianalytic expression for the piston 
removed phase covariance on a circular aperture in the presence of Komolgorov turbulence (Tyler 
1994). The covariance between the piston removed wavefront phase from two different sources at 
two different points fi and rb in the telescope pupil plane is given by 
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Here D is the aperture diameter, k is the wavcnumber, and C^{z) is the turbulence profile along 
the line of sight to the star. The numerical constant S is given by 
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where F (z) is the Gamma function. The quantity Qab (z) is 

^abiz) = (-^j^ab 



(8) 



where 9ab is the angular offset between the two sources. This parameter may be interpreted as the 
altitude dependent shear between the two beams. Finally, the two functions Gi and G2 in Equation 
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6 are defined in terms of the Gauss hypergeometric function 

T{a + n)T{b + n)T{c) 
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Note that the covariance expression contains dependencies on observing wavelength, aperture di- 
ameter, turbulence profile, and angular offset between the two stars. The range variable z implicitly 
incorporates the dependence of the covariance on zenith angle, which enters through C^{z), ^ab{z), 
and the range integral itself. 

The expression for the covariance in Equation 6 is of great utility, and a number of familiar 
results in adaptive optics may be derived from this expression. Three limiting cases are shown 
in Appendix A: the phase structure function in the presence of uncompensated turbulence, the 
aperture averaged phase variance in the presence of uncompensated turbulence, and the aperture 
averaged residual phase variance due to anisoplanatism. These cases motivate the definition of the 
Pried parameter ro and the isoplanatic angle which are restated in this Appendix. 

Expanding the expression for the anisoplanatic structure function using Equation 2 permits 
expression of -Dapi(n_) ^2) as a sum over ten phase covariance terms involving (pai^ aiid (/)fo(r). Direct 
evaluation of the resulting expression using Equation 6 yields 
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This result indicates that the anisoplanatic structure function Dapi (rt, ^^2) is in fact only a function 

of fi — r2, and is therefore stationary over the pupil plane. This property permits the factorization 
of the anisoplanatic and guide star OTFs shown in Equation 5. 

As in the case of the covariance, I?apl (^) depends on observing wavelength, aperture diameter, 
turbulence profile, zenith angle and angular offset between the guide star and the direction of 
interest. Examples of the anisoplanatic structure function are shown in Figure 1 for a particular 
set of these parameters. 
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Fig. 1. — Anisoplanatic structure functions for a star at a 20 arcsec offset from the guide star, 
computed using Equation 12. The left and right panels show Dapi(^; 0) at observing wavelengths of 
2.12 fim and 1.65 /xm as a function of location r within the pupil plane of a 5 meter telescope. The 
angular offset from the guide star was oriented along the Y axis, and the broken X-Y symmetry 
gives rise to the anisotropy apparent in these functions. This anisotropy is responsible for the 
radial elongation of the off axis PSF that is a characteristic feature of anisoplanatism. Despite this 
anisotropy the structure function is stationary over the pupil plane, so that -Dapi(n_5 "^2) = -Dapi(r2 — 
n). Typical atmospheric turbulence conditions at Palomar were assumed in these calculations. 

2.3. The Anisoplanatic OTF and PSF 

Equation 5 represents a model for the OTF at any point in the field of view. This quantity is a 
strong function of many parameters, which enter through both the anisoplanatic structure function 
Da,pi (r) and the guide star OTF. These dependencies give rise to the rich morphology in the 
adaptive optics PSF that presents such a challenge in the analysis of adaptive optics observations. 

The classic dependence displayed by the adaptive optics PSF is its degradation with increasing 
angular offset from the guide star. Figure 2 contains a series of plots that illustrate this behavior 
for a 5 meter telescope under typical atmospheric turbulence conditions at Palomar. The first 
plot shows the aperture averaged residual RMS optical path difference (OPD) vs angular offset 
from the guide star. These results were computed by evaluating ([4>apii^)]^^ using Equation 6, 
and averaging this quantity over the pupil plane. The OPD is found by computing the square 
root of this result and dividing by k. Also shown in this plot is the OPD computed from the 
approximation described in Appendix A. Even at small angular offsets, significant discrepancies 
exist between the exact result and the approximation. 

Consider an idealized adaptive optics system that acts to eliminate all wavefront error in the 
direction of the guide star, so that (j){f) = 0. In this case the OTF in Equation 5 depends only 
on the anisoplanatic structure function -Dapi ("f^, which may be evaluated using Equation 12. This 
allows computation of an OTF that incorporates only the effects of anisoplanatism, from which a 
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Fig. 2. — Illustrative effects of ansioplanatism on a 5 meter telescope, as described in the text. The 
upper left panel shows the aperture averaged residual RMS OPD due to anisoplanatism vs angular 
offset. The upper right panel shows the 2.12 Strehl ratio vs angular offset. The lower panels 
show the 2.12 //m anisoplanatic OTF and PSF at an angular offset of 20 arcsecs from the guide 
star. The PSF has been normalized by the peak value of the ideal, diffraction limited PSF. Typical 
atmospheric turbulence conditions for Palomar were assumed for these calculations. 

PSF may be evaluated by Fourier transformation. Figure 2 shows the 2.12 /im anisoplanatic OTF 
and PSF at an angular offset of 20 arcsecs from the guide star, computed using this method. The 
PSF displays an increase in elongation along the direction to the guide star that is a characteristic 
of anisoplanatism. At this angular offset, the Strehl ratio has dropped to 47% due solely to the 
effects of anisoplanatism. 

In the same way, the anisoplanatic Strehl ratio may be computed as a function of angular offset 
from the guide star. These results are also plotted in Figure 2. For comparison, the anisoplanatic 
Strehl ratios computed from the Marechal approximation are also shown. In this approximation, 
the aperture averaged phase variance computed from the approximation was exponentiated to 
form an estimate of the anisoplanatic Strehl ratio. The combination of these two approximations 
significantly underestimates the anisoplanatic Strehl ratio computed directly from the anisoplanatic 
PSF. 
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This approach to evaluating the anisoplanatic OTF and PSF is illustrative, but the application 
of Equation 5 to observational data requires a treatment of the guide star OTF. While significant 
efforts have been made to estimate the guide star OTF using the statistical properties of 0a (^) 
(Veran et al. 1997), evaluation of Equation 5 is considerably simplified if an observation of the 
guide star may be used in the prediction of the PSF elsewhere in the field of view. For this pro- 
cedure to work, the guide star must itself be a point source. This is often the case in adaptive 
optics observations. In these circumstances, the guide star PSF may be extracted from the ob- 
servational data and Fourier transformed to form the guide star OTF that appears in Equation 
5. The anisoplanatic structure function Dapl (0 may be evaluated from Equation 12 for the field 
point of interest. The only parameter in Equation 12 that is not defined by the observation is the 
turbulence profile, which may be measured using an independent set of equipment. The product of 
the observed guide star OTF and the anisoplanatic transfer function provides a prediction of the 
OTF at the field point of interest. Fourier transformation of this OTF yields a prediction for the 
PSF. There are no free parameters in this prediction. 



3. Observations 

To carry out a comparison of the model in Equation 5 with AO compensated image data, a set 
of observations was carried out at Palomar Observatory on the night of August 20, 2005. In this 
experiment, a well separated binary system was chosen and one member of the binary was used 
as the guide star for the Palomar Adaptive Optics system (Troy et al. 2000) on the Hale 5 meter 
telescope. A sequence of AO compensated near infrared images of the binary system from the 
PHARO infrared camera (Hayward et al. 2001) and a sequence of turbulence profiles from a set of 
DIMM/MASS equipment were acquired contemporaneously. Using the methodology presented in 
Section 2, a prediction for the PSF of the binary companion was formulated from the observed PSF 
of the guide star, the measured turbulence profiles, and the parameters defined by the observation. 
This prediction was then compared directly to the observed images of the binary companion. This 
section describes the details of these observations. 



3.1. Turbulence Profile Measurements 

Turbulence profile measurements were acquired using a DIMM/MASS unit on loan from the 
Thirty Meter Telescope Project. The unit itself has been described by Skidmore et al. (2004), 
and consists of a 35 cm robotic telescope that feeds both a DIMM and MASS instrument. The 
methodology employed in deriving turbulence profiles from the combination of DIMM and MASS 
measurements is discussed by Tokovinin et al. (2003) and references therein. The DIMM measures 
the differential motion between images acquired through two adjacent apertures. This measurement 
is sensitive to the integrated turbulence profile. A MASS unit measures the scintillation index of 
images acquired through four apertures of different radii. These measurements are sensitive to the 
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Fig. 3. — Turbulence parameters on the night of August 20, 2005 measured using the DIMM/MASS 
equipment at Palomar. These parameters were computed from the seven layer profiles estimated 
from the DIMM/MASS measurements, and are quoted at zenith and at a reference wavelength of 
.5 ^m. The mean value of the Fried parameter tq was 7.1 cm, and that of the isoplanatic angle 
^0 was 2.54 arcsecs. Both parameters displayed considerable excursions from their means over the 
course of the night. 

distribution of higher altitude atmospheric turbulence, which gives rise to scintillation. Together 
these measurements permit an estimate of the turbulence profile at altitudes of 0, .5, 1, 2, 4, 8 and 
16 km. In a comparison of the turbulence strength measured at each of these layers, simultaneous 
MASS and scidar observations on Mauna Kea have displayed agreement to a factor of several 
(Tokovinin et al. 2005). One of the difficulties in the reconstruction of the turbulence profile from 
MASS data is that the procedure is somewhat ill-conditioned, in that simultaneously adjusting the 
layer altitude and the strength of turbulence can lead to similar scintillation indices. However, to 
leading order anisoplanatism displays this same type of covariance (c.f. Equation A13). In this 
sense both anisoplanatism and MASS are low resolution atmospheric profilometers. 

The DIMM/MASS at Palomar Observatory was located 300 meters from the Hale 5m telescope. 
This unit was installed on the roof of a building approximately 10m above ground level. Turbulence 
profiles were measured about once every 90 seconds throughout the night of August 20th. For the 
first hour of the experiment, the DIMM/MASS unit tracked Beta Drac at 17''30"' + 52°18". A 
half hour gap in coverage occurred between 5:00 and 5:30 UT due to a tracking error experienced 
by the robotic telescope. Prom 5:30 UT onwards, the unit recorded images from Alpha Cep at 
2l'*19"* + 62°35". For comparison, the binary imaged using the AO system was at 17'*59"^ + 64°08". 

Figure 3 shows the values of ro and computed from these profiles. Both of these parameters 
display deviations from their means of order a factor of two over the course of the night. The 
variability in the turbulence profiles imply a factor of two variation in the RMS OPD due to aniso- 
planatism. This in turn yielded pronounced variability in the Strehl ratio of the binary companion 
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that was readily detectable in the AO compensated image data. 

3.2. Adaptive Optics Observations 

Adaptive optics observations of the binary system HD164983 + HD164984 were acquired over 
a three hour period between 4:17 and 7:18 UT. This binary system has an angular separation of 
21" and is oriented at 282° east of north. The Johnson V magnitudes of HD164983 and HD164984 
are 7.6 and 7.3, respectively, and the latter source was used as the guide star for the adaptive 
optics system. The 25 arcsec field of view of the PHARO camera was used, permitting both 
binary members to be positioned on the detector simultaneously. Imaging was performed using 
two narrowband filters: an H2 filter with a central wavelength of 2.123 fim and a bandpass of .007 
fj,m and an Fell filter with a central wavelength of 1.648 /im and a bandpass of .03 ^m. Exposure 
times of 2.8 seconds and 1.4 seconds were chosen for the H2 and Fell observations, respectively. 
The resulting stellar images peaked at less than 30% of the detector full well depth, which is well 
within the linear range of the detector. Images were acquired in a 7 point linear dither pattern, in 
which the binary was shifted up and down the detector in 4.5" steps. Ten exposures were acquired 
at each dither position. After culling the data for bad images, a total of 703 exposures at 2.123 /xm 
and 384 exposures at 1.648 /xm remained. 

Calibration of the image data was carried out in the customary way. Flat field calibration was 
performed using twilight sky flats. Sky subtraction was performed by forming the median of the 
dithered exposures and subtracting this median from each exposure. Finally, 3.8" subimages cen- 
tered on each of the two binary members were extracted from each exposure for use in the analysis 
below. The angular extent of these subimages was chosen to encompass all residual scattered light. 

4. Data Analysis 

A prediction of the PSF at the location of the binary companion was formulated for each 
exposure using the method described in Section 2. Each subimage of the guide star was Fourier 
transformed to form the guide star GTE. An anisoplanatic GTE was computed for each measured 
turbulence profile using the observing wavelength, aperture diameter, angular offset to the binary 
companion, and zenith angle at the time of the measurement. For each adaptive optics exposure, 
an anisoplanatic GTE was formed for the time of the exposure by interpolating the anisoplanatic 
OTEs computed for the nearest turbulence profile measurements. The pointwise product of the 
guide star OTF and the interpolated anisoplanatic OTE was formed, and the resulting GTE was 
Fourier transformed to generate a PSF prediction for the binary companion sampled at the 25 mas 
pixel scale of the PHARG image data. 

For each exposure, the Strchl ratios of each observed binary member and of the predicted 
companion PSF were calculated. To compute these Strehl ratios, an ideal, diffraction limited PSF 
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for the Hale 5m was simulated by forming the pupil plane wavefront, Fourier transforming to form 
the image plane wavefront, and computing the square modulus. The pupil plane wavefront used 
in this procedure accounted for the shadows cast by the secondary mirror and the four struts 
that support it. For each subimage this diffraction limited PSF was normalized to have the same 
integrated signal, and the Strehl ratio was computed as the ratio of the peak value in the image 
to that of the ideal PSF. In these observations the 2.12 ij,m images were oversampled at 1.7 times 
Nyquist, while the 1.65 iim images were oversampled at 1.3 times Nyquist. This oversampling 
mitigates many of the subtle effects that can arise in computing Strehl ratios from image data 
(Roberts et al. 2004). 

In addition to the Strehl ratio analysis, a direct comparison was carried out between the 
predicted and observed PSFs of the binary companion. Because the guide star PSF was used in 
formulating this prediction, the predicted PSF represents the image that would be obtained for 
a point source with the same brightness as the guide star. In fact the binary companion has a 
different brightness, which may be varying in time. Possible sources of such variability include the 
presence of cirrus clouds during observations or intrinsic photometric variability of the star itself. 
In addition, subimages were extracted from the PHARO cxposTircs at a specific angular offset. 
Variations in the binary offset may occur from one exposure to the next, and it is of interest to 
measure these offsets from the data. These parameters constitute the differential photometry and 
astrometry of the binary, and their measurement requires a fit of the predicted companion PSF to 
the observed companion PSF. 

For this analysis, a simple four parameter model was used. The relationship between an 
observation of the binary companion P{x,y) and the predicted PSF R{x,y) was taken to be 

P{x,y) = bR{x - Ax,y - Ay) + c (13) 

Here b represents the differential amplitude between the two binary members, while Ax and Ay 
represent the differential angular offsets. The constant c models any differential in the background 
level of the guide star and companion data, and guarantees that the residuals from the fit have zero 
mean. Solution of this equation is readily performed using the predicted and observed companion 
OTFs. In OTF space, the parameters Ax and Ay induce a phase slope by the shift theorem 
(Bracewell 1986). A merit function is readily formulated in this space, and minimization of this 
function leads to a simple iterative solution for all four parameters. 

The outcome of this fitting procedure permits measurement of the differential astrometry 
and photometry between binary members. The overall angular offsets assumed in the subimage 
extraction were added to the fitted values of Ax and Ay to yield the differential astromctric offset 
between the binary members. The differential photometry was computed as the ratio of the total 
flux of the fitted PSF to that of the guide star. Because the residuals are zero mean, the total 
flux of the fitted PSF is identical to that of the observed image of the binary companion. At first 
glance, the use of the fitted PSF appears to provide no advantage over the observed PSF. However, 
these images were oversampled and their OTFs were lowpass filtered to suppress spatial frequencies 
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Fig. 4. — Observations of the guide star HD164984 and its companion HD164983 at 2.12 fim 
and 1.65 fim. These observations were formed by integrating five sequential exposures, and are 
displayed on a log stretch. The angular offset between these stars is 21.3 arcsecs, and only a 2 
arcsec field surrounding each star is shown. The first and second columns show observations of the 
guide star and companion, respectively. The predicted PSF of the companion appears in the third 
column, and was formulated from the product of the guide star OTF and the anisoplanatic OTF 
as described in Section 2. The Strehl ratios are shown at the bottom of each image. In these data, 
anisoplanatism has degraded the Strehl ratio of the companion by factors of two to three relative 
to the guide star. This degradation is accurately captured in the predicted companion PSF. 

larger than the cutoff set by the telescope aperture. This filtering step substantially reduces the 
noise in the image. 

5. Results 

Examples of the 2.12 fj,m and 1.65 /um observations of the binary members are shown in Figure 
4. These observations were integrated over five sequential exposures, and serve to illustrate the 
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quality of the AO compensation delivered by the PALAO system. In the particular 14 second 
observation at 2.12 ^um shown in this figure, the guide star Strchl ratio was 51%. Due to the effects 
of anisoplanatism, the measured Strehl ratio of the binary companion was degraded to 28%. The 
predicted companion PSF computed from the guide star OTF and the turbulence profile at the 
time of this observation is also shown in this figure. The Strehl ratio calculated for the predicted 
companion PSF was 27%. The level of agreement between measured and predicted companion 
Strehl is a strong indication that the OTF model in Equation 5 accurately captiircs the effects 
of anisoplanatism. This also implies that the turbulence profiles measured by the DIMM/MASS 
equipment and used in this model reflect the true distribution of atmospheric turbulence. A similar 
level of agreement is seen for the 7 second observation at 1.65 fim shown in the same figure. 

Figure 5 shows the values of tq and computed from the measured turbulence profiles that 
were acquired over the course of the three hour observation. These parameters have been computed 
for the zenith angle of the guide star, which varied between 30 and 40 degrees during the observa- 
tions. The values of these turbulence parameters display considerable temporal variation. Figure 5 
also shows the time dependence of the Strehl ratios at 2.12 //m. These Strehl ratios have again been 
averaged over five sequential exposures. The measured guide star Strehls vary significantly, and 
show little correlation with the Fried parameter. This suggests that sources of wavefront error other 
than residual fitting error contributed significantly to the guide star error budget. The measured 
Strehl ratios of the binary companion are also plotted, along with the Strehl ratios computed from 
the predicted companion OTF. As in Figure 4, these predictions are in excellent agreement with the 
measurements, and arc able to track the companion Strchl ratio to a few percent despite a factor 
of two variability in the measured Strchl ratio of both the guide star and the companion. Also 
plotted are the Strehl ratios of the binary companion predicted using a traditional error budget 
approach. These predictions were generated from the product of the guide star Strehl ratio and the 
anisoplanatic Strehl ratio computed using the approximation, as described in Section 2.3. These 
predictions fall well below the measurements, as expected from the discrepancies in the anisopla- 
natic Strehl ratio shown in Figure 2. Lastly, Figure 5 contains a plot of measured vs. predicted 
2.12 /im and 1.65 fim Strchl ratios for the binary companion. This plot again illustrates strong 
agreement between the measured Strehl ratios and those derived using the predicted companion 
OTF. 

Results from fitting the predicted companion PSF to the data are shown in Figure 6 for the 
observations shown in Figure 4. This figure shows that the prediction of the companion PSF slightly 
overestimates the core and underestimates the wings of the observed PSF. Horizontal cuts through 
the simulated diffraction limited PSF display complex morphology arising from light scattered by 
the four struts supporting the secondary mirror. The scattered light is most pronounced along the 
horizontal and vertical axes of the image, which are aligned with the support struts. Horizontal 
cuts through the residuals indicate that the predicted companion PSF agrees with the observations 
to about 10% accuracy at any point in the image. Stated another way, this PSF fitting procedure 
has improved the dynamic range of the observation by an order of magnitude. 
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Fig. 5. — Turbulence parameters and Strehl ratios over the course of the observations. The upper 
panels show the Fried parameter and isoplanatic angle at .5 computed for the zenith angle of 
the guide star, which ranged from 30 to 40 degrees over the three hour observation. The lower left 
panel shows the measured 2.12 //m Strehl ratios of the guide star and companion. Also plotted 
are the Strehl ratios predicted from the OTF formulation in Section 2, and predicted using the 
^0 approximation described in the text. The lower right panel shows the predicted vs. measured 
Strehl ratio of the binary companion. Strehl ratios at both 2.12 ^m and 1.65 ^m are included 
in this plot. These plots indicate that the OTF formulation accurately predicts the Strehl ratio 
degradation that arises from ansioplanatism. See the electronic edition of the Journal for a color 
version of this figure. 

These residuals display two regimes, in which the accuracy of the predictions was limited by 
different effects. At radii less than about .5 arcsecs, the residuals were dominated by systematic 
errors between the model and the observed data. These systematics may arise from a number of 
different effects. The power spectrum of atmospheric turbulence may be non-Komolgorov, so that 
Equation 12 is only approximately correct. Accuracy in the turbulence profile measurements from 
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Fig. 6. — Results of fitting the model of the companion PSF to the observed data for the 2.12 /im 
and 1.65 fim observations shown in Figure 4. The left column shows the residuals in a 2 arcsec 
field around the companion, displayed on a log stretch. The right column shows horizontal cuts 
through the simulated diffraction limited PSF, the observed guide star PSF, and the predicted and 
observed companion PSF. The complex behavior displayed in the wings of the diffraction limited 
PSF arises from light scattered by the struts supporting the secondary mirror. The magnitude of 
the residual difference between the predicted and measured companion PSF is also plotted. At any 
point within the image, the predicted companion PSF matches the observed data to an accuracy 
of about 10 percent. See the electronic edition of the Journal for a color version of this figure. 
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the DIMM/MASS unit may limit the quality of the prediction. Finally, nonlinear response or charge 
diffusion in the infrared detector in PHARO may generate differential errors between the measured 
guide star and companion PSFs. At an angular separation of about .5 arcsccs, the signal level in the 
residuals dropped below a noise floor. This noise floor is set by a combination of the detector read 
noise, shot noise and quantization noise in the analog to digital conversion. The dominant source 
of noise depends on the efficacy of the lowpass spatial filtering that was performed on the OTF, 
which itself depends on the degree to which the images were oversamplcd. Further investigation 
will be required to determine the noise source that limits the precision of this technique in these 
two regimes. 

At a radius of .5 arcsecs the value of the residual intensity was 3 to 4 orders of magnitude 
less than the peak of the measured companion. This level of PSF rejection is comparable to that 
achieved in near infrared AO observations that employ a Lyot coronograph. These systems provide 
about 4 to 5 orders of magnitude rejection at an offset of 1" from a star (Oppenheimer et al. 2000; 
Nakajima et al. 2005). In comparing Lyot coronography to this fitting technique, a disadvantage 
of the latter is that it requires the use of another star to serve as the PSF reference. On the other 
hand, Lyot coronographs employ an opaque focal plane mask that limits the inner working radius 
of the observation. For example, the Lyot coronograph in PHARO has masks with diameters of 
.41 and .91 arcsecs. There is no inner working radius in this fltting technique. 

Figure 7 shows the differential photometry between the binary members derived from fitting 
the predicted companion PSF to the observed data. The timeseries of differential photometric 
measurements shows periods of substantial variability at 05:30 UT for the 2.12 /xm observations, 
and at 07:15 UT for the 1.65 /xm observations. The origin of this variability is not clear, but could 
plausibly be ascribed to cirrus. When these data were excluded, the standard deviation of the 
measured differential photometry was 1 part in 10^ for both the 2.12 fim and 1.65 fim exposures. 
The mean differential photometry between the binary members is shown in column 2 of Table 1. 
The uncertainties quoted on these measurements are the estimated errors of the mean (Bevington 
& Robinson 1992), and are less than 1 part in 10^. 

This level of photometric stability may be compared to recent results. Roberts et al. (2004) 
obtained differential photometric stability of order 1 part in 10^ in J, H, and K band adaptive 

Table 1. Differential photometry and astrometry for the binary system HD164984+HD164983. 



A Differential a\\ a± 
(/xm) Photometry (arcsec) (arcsec) 



P 

(arcsec) 



P.A. 

(deg) 



2.12 .7903±.0004 .0152 .0100 21.3322±.0006 282.3923±.0011 
1.65 .8316±.0007 .0228 .0154 21.3306±.0012 282.3932±.0023 
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Fig. 7. — Differential photometric and astrometric measurements of the binary 
HD164984+HD164983 for 2.12 fxm (upper) and 1.65 /xm (lower) observations. Each point 
in these plots corresponds to a single 2.8 second exposure at 2.12 /xm, or 1.4 second exposure at 
1.65 /im. The left panels show the measured flux ratio of the companion to primary throughout 
the course of the observations. The right panels show the measured differential angular offset 
between the binary members for each exposure in the experiment. The elliptical scatter in the 
astrometric data arises from differential atmospheric tilt jitter between the two stars that is 
induced by atmospheric turbulence. This tilt jitter is predicted to be larger along the orientation 
of the binary, which lies along the xy axis indicated in these plots. 

optics observations of several binary systems on the Advanced Electro-Optical System (AEOS) 3.6 
m telescope. Seeing limited observations in the optical at the 6.5 m Multiple Mirror Telescope 
(MMT) (Hartman et al. 2005) and in the near infrared at the 3.8 m United Kingdom Infra-Red 
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Telescope (UKIRT) (Snellen 2005) displayed differential photometric stability of 1 part in 10^. 
All of these observations employed broader filters and lower total integration times than were 
used in this experiment. Relative to seeing limited observations, the sensitivity gains provided by 
adaptive optics have afforded this level of photometric stability at much lower flux levels, allowing 
observations of fainter targets. This technique may be usefully applied to observations of eclipsing 
binaries, transiting planets, and other systems that display photometric variability in the near 
infrared. 

Figure 7 also shows the differential astrometry derived from the fit. Each exposure yielded 
a single measurement, aU of which were combined to form the scatter plots in this figure. The 
astrometric offsets display an elliptical scatter, with larger errors along the axis connecting the 
binary members. This behavior is consistent with that of differential atmospheric tilt jitter, which 
arises from anisoplanatism of the tilt component of atmospheric turbulence. This effect leads to a 
random, achromatic fluctuation in the relative displacement of two objects. The standard deviation 
of this tilt jitter differs along the axes parallel and perpendicular to the orientation of the binary. 
A three term approximation to the parallel and perpendicular components of the variance arising 
from differential atmospheric tilt jitter is given by Sasiela (1994). 
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Here the turbulence moments /i^ are defined as 
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(14) 



(15) 



Using the mean turbulence profile over the three hour observation, the parallel and perpendicular 
components of differential atmospheric tilt jitter were computed from Equation 14 to be a\\ = 36 
mas and ctj_ = 21 mas. These values are somewhat larger than those measured from the astrometric 
data, which appear in columns 3 and 4 of Table 1. The discrepancy likely arises from the finite 
integration time of the exposures. Assuming a characteristic wind speed of 5 m/s, these integration 
times are comparable to the wind crossing time for the 5 meter aperture, suggesting that tilt jitter 
has partially averaged away. This hypothesis is further supported by the reduction in tilt jitter 
between the 1.4 second exposures at 1.65 /Ltm and the 2.8 second exposures at 2.12 /Ltm. 

The mean differential astrometry between the binary members is shown in columns 5 and 
6 of Table 1. A PHARO pixel scale of .02522 asec/pixel was used in the calculation (Metchev 
& Hillenbrand 2004). The uncertainties quoted for these measurements are again the estimated 
errors of the mean. This uncertainty is about 1 mas in the binary separation p, and a few arcsecs 
in position angle. Note that the differential astrometry of the 2.12 /xm and 1.65 jim observations 
are in agreement at about the la level. This provides an independent validation of the accuracy of 
these measurements. 



This level of accuracy may be compared to recent astrometric results obtained with the adaptive 
optics system on AEOS. Roberts et al. (2004) performed I band adaptive optics observations of a 
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number of binaries and measured the differential astrometry to accuracies of 10 to 20 mas. These 
binaries had separations up to 5 arcsecs. Drummond ct al. (2003) reported H band observations 
of a 400 mas binary with an astrometric precision of 1 mas. Note that to leading order the tilt 
variance in Equation 14 scales as 6"^ /D'^/^, so that these binaries suffer much less from differential 
atmospheric tilt jitter than the 21 arcsec binary observed in this experiment. Compared to the 
observations at Palomar reported here, both of these experiments employed broader filters and 
lower total integration times. 

These results may also be compared to those obtained using other astrometric techniques. An 
astrometric accuracy of order 100 jias has been achieved on the Palomar Testbed Interferometer 
(PTI) for a 30 arcsec binary (Lane et al. 2000), and recently accuracies of tens of //as have been 
obtained for binaries with separations less than an arcsec (Lane & Muterspaugh 2004; Muterspaugh 
et al. 2006). This instrument has a limiting magnitude of about 6 due to the 40 cm size of its 
apertures. The STEPS program on the Hale 5m telescope at Palomar has achieved astrometric 
accuracies of less than 1 mas using seeing limited observations in the visible. (Pravdo & Shaklan 
1996). These observations employed a 2 arcminute field of view and imaged crowded fields so as 
to establish an astrometric reference grid from the multiple objects in the field. Both of these 
techniques employ much broader filters than were used in the observations reported here. Further 
experiments will be required to understand the circumstances under which the astrometric technique 
described in Section 4 is competitive with these other methods. Generally speaking, the sensitivity 
improvements afforded by adaptive optics on a 5 meter aperture will permit application of this 
methodology to fainter limiting magnitudes than are accessible with PTI or STEPS, substantially 
increasing the number of accessible targets. Differential atmospheric tilt jitter scales as D~'^/^ , and 
application of this technique on a larger aperture telescope will reduce this jitter while at the same 
time increasing the signal to noise ratio and decreasing the width of the PSF core. This would 
provide substantial improvements on the astrometric accuracies that have been reported in this 
experiment. 



6. Conclusions 

The research presented here has drawn together a number of different elements in order to 
generate predictions of the adaptive optics PSF. Factorization of the OTF in Equation 5 permits 
one to use the guide star PSF as a reference for observations throughout the field of view. This 
PSF encapsulates the complex behavior of the adaptive optics system that is otherwise very dif- 
ficult to model. The covariance expression in Equation 6 provides an analytic formulation that 
captures the dependencies of anisoplanatism on aperture diameter, observing wavelength, turbu- 
lence profile, angular offset and zenith angle. Measurements of the turbulence profiles from the 
DIMM/MASS equipment provide the one input parameter for these predictions that is not deter- 
mined directly from the observations. These three elements provide a methodology for adaptive 
optics PSF prediction that is accessible to direct experimental validation. 
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The binary star observations described above are in excellent agreement with these predictions. 
The Strehl ratios computed from the predicted companion PSFs match the measured values to an 
accuracy of a few percent, despite a factor of two temporal variability in both the guide star and 
companion Strehl ratios. The predicted companion PSF matches observations to about 10% out 
to radii of 1 arcsec. While this agreement serves to validate the predictive methodology, it must 
be emphasized that the adaptive optics PSF depends on a large number of parameters. This 
three hour experiment has tested these predictions over a very modest region of the underlying 
parameter space. A broader application of this methodology at shorter observing wavelengths, over 
wider fields of view, and under diverse turbulence conditions will provide a more stringent test. In 
some observing conditions these predictions will almost certainly fail due to the approximations 
discussed in Section 2. Likewise, the target of this experiment was a relatively bright binary whose 
members are of nearly identical magnitude. Additional observations will be required to understand 
the degree to which the photometric and astrometric precision and the contrast levels reported 
here are attainable in more diverse observational programs. Experiments at Palomar and Keck 
Observatories arc currently being planned to perform this experiment on more binaries, and to 
employ this methodolgy on crowded field image data using both natural and laser guide star AO 
systems. 

Agreement between the predicted and observed companion PSF indicate a level of consistency 
between turbulence profile measurements from the DIMM/MASS equipment and the effects of 
anisoplanatism on the AO compensated image quality. Both Multiaperture Scintillation Spectrom- 
etry and anisoplanatism are sensitive to higher altitude turbulence, but each is sensitive to the 
underlying turbulence statistics in a different way. The former is sensitive to aberrations intro- 
duced by short spatial wavelengths, which give rise to scintillation. In contrast, anisoplanatism 
arises from aberrations at all spatial wavelengths, though the manner in which these aberrations 
contribute depends on the altitude of the turbulence and the spatial frequency of the aberration. 
In this way, these two types of measurements are sampling different spatial frequency ranges of the 
turbulence power spectrum. At some level the consistency of the predictions and observations serve 
to validate the assumption of a Komolgorov turbulence spectrum. Further efforts will be required 
to understand the level of agreement implied by these and future results. One possible approach 
is to perform a sensitivity analysis by generalizing Equation 6 to include other classes of power 
spectra. An analysis of this type has been carried out by Lazorenko (2002), who considered the 
effects of non-Komolgorov power spectra on differential atmospheric tilt jitter. 

These results suggest several lines of longer term development that may be of direct benefit to a 
number of astronomical applications. The plot of Strehl ratio vs. guide star offset shown in Figure 2 
illustrates that the benefits of adaptive compensation occur over fields much larger than anticipated 
from the approximation. This suggests that adaptive optics systems with fields of view of order 
several arcminutcs may be usefully employed for near infrared observations. One of the challenges 
in interpreting these data over wide fields arises from the temporal and field dependent evolution 
of the adaptive optics PSF. The methodology described above allows a quantitative analysis of 
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such wide field observations that can account for these effects. The photometric and astrometric 
results presented in Section 5 and the dynamic range improvements implied by Figure 6 serve as 
illustrations of the astronomical potential of this methodology. Application of this technique to 
the imaging and deconvolution of crowded fields and extended objects would constitute a natural 
progression. The approach could readily be applied to the deconvolution of both image data and 
spatially resolved spectra acquired with an integral field unit. 

Lastly, agreement between predictions and observations serve as an important on-sky vali- 
dation that anisoplanatism is accurately understood in the context of near infrared astronomical 
observations over arcminute fields of view. Anisoplanatism constitutes the fundamental process 
underlying the use of tomography in adaptive optics. Tyler (1994) describes the application of the 
phase covariance in Equation 6 in combining wavefront measurements from multiple guide stars 
to form an estimate of the wavefront in a different direction. That this same expression has been 
used in this research to accurately predict the anisoplanatic degradation of the near infrared adap- 
tive optics PSF over an arcminute field is a strong indication that tomographic algorithms will be 
successful in the same observational context. 
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A. Three Familiar Results in Adaptive Optics 

This Appendix uses Equation 6 to recover three familiar results in adaptive optics. These 
results serve to illustrate the validity of Equation 6, and its broad applicability in performing 
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analytic and numerical calculations in adaptive optics. Equation 6 represents the piston removed 
phase covariance on a circular aperture in the presence of Komolgorov turbulence, and the three 
results below are valid under these assumptions. 



A.l. Phase Structure Function for Uncompensated Turbulence 

As a first example, consider the phase structure function (ri, r2) in the presence of uncom- 
pensated turbulence. 

D^in,r2) = {[<Piri)-4>ir2)f) (Al) 
= ([<!> {ri)f) + ([^ {r2)f) - 2 (0 (ri) (r^)) 



The three covariance functions may be rewritten in terms of Equation 6. In the resulting expression, 
all dependencies on the functions Gi and G2 drop out in the difference, leaving only 

D^iri,r2) = 2'/'Ek^\{n-f2)\'/^ JdzCliz) (A2) 

Define the Pried parameter ro as 

r-'/' = 2'/'^eJ^ dzCliz) (A3) 



where the constant A is 

A 



-| 5/6 

= 6.88388 (A4) 



The phase structure function may be rewritten in terms of ro as 

i^^(n,f2) = A(^^i-^J (A5) 

This is the well known expression for the piston removed phase structure function on a circular 
aperture in the presence of Komolgorov turbulence. 



A. 2. Aperture Averaged Phase Variance for Uncompensated Turbulence 



As a second example, the aperture averaged phase variance in the presence of uncompensated 
turbulence is evaluated. This quantity may be computed by integrating the phase variance <^['/'(r)]^^ 
over the circular aperture and dividing by the area of the aperture. The phase variance may be 
rewritten using Equation 6, yielding 
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The integral over the hypergeometric function in Gi may be performed term by term, yielding 
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where the second equality has employed the relationship (Abramowitz & Stegun 1972) 

r(c)r(c -a-b) 



2F1 (a, 6; c; 1) = 

The second term in Equation A6 may be written 
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The integral has been evaluated analytically using the substitution y = cos^. Further manipula- 
tions involving relationships between Gamma functions or direct numerical calculation show that 
the first term evaluates to exactly twice the second. The resulting expression becomes 
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This is the well known result for the aperture averaged phase variance in the presence of uncom- 
pensated turbulence. 



A. 3. Aperture Averaged Residual Phase Variance from Anisoplanatism 



As a final example of particular relevance to this paper, consider the aperture averaged residual 
phase variance due to anisoplanatism. This quantity may be computed by integrating the residual 
phase variance ^[0fo(r) — 0a(r)]^\ over the circular aperture and dividing by the area of the aperture. 
Again using Equation 6, we find 
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The first term is 
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where the isoplanatic angle is defined as 

^-5/3 ^ 28/3sA:2 / dz Cl{z)z'/^ 
Jo 



(A13) 



This is the well known isoplanatic angle approximation to the aperture averaged residual phase 
variance due to anisoplanatism. 

Consider expanding the remaining three pairs of terms in Equation All in a Taylor series about 
^ab{^) = 0- The constant terms cancel in the differences. The linear term in Clabi^) vanishes in the 
function G2 ( ^ab{z) ) , since it is a function only of the magnitude of this vector. The linear term in 

since 



^ab{^) also vanishes in the sum of the terms Gi ^ 2r/-D — ^ab{z) ^ and Gi ^ 2f/D + ^ab{z) ^ , 

this vector enters these two terms with opposite sign. Thus, the remaining terms in Equation All 

2 



have a leading order dependence of 



These terms are only slightly higher order than the 
first term in this equation, and their contribution to the aperture averaged residual phase variance 
can be significant even at modest angular offsets. This is illustrated by the large discrepancy 
between the isoplanatic angle approximation to the aperture averaged residual phase variance and 
the exact result shown in Figure 2. 
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